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Abstract. - We study numerically how conservation laws affect the optical conductivity 
(t(ij;) of a slightly doped one-dimensional Mott insulator. We investigate a regime where the 
average distance between charge excitations is large compared to their thermal de Broglie wave 
length and a classical description is possible. Due to conservation laws, the dc-conductivity is 
infinite and the Drude weight D is finite even at finite temperatures. Our immerical results 
test and confirm exact theoretical predictions for D both for integrable and non-integrable 
models. Small deviations from integrability induce slowly decaying modes and, consequently, 
low-frequency peaks in o-((jj), which can be described by a memory matrix approach. 



Conservation laws and slowly decaying modes strongly affect transport properties . The 
reason is that the presence of conserved quantities (e.g. momentum) can prohibit the complete 
decay of an electrical current J (in the absence of external fields). If a certain fraction 
Jc of J does not decay, the dc-conductivity is infinite and the optical conductivity a(uj) is 
characterized by a finite Drude weight D on top of a regular contribution ai-og{uj) even at 
finite temperature T, Iiea{Lu) — 2TrDS{uj) + Reo'rcgl'^)- 

In quasi one-dimensional (ID) systems conservation laws are especially important. Even 
in the presence of Umklapp scattering from a periodic potential, certain pseudo-momenta are 
approximately conserved Further, the low-energy properties of interacting electrons in 
one dimension are well described by integrable models like the Luttinger, Sine-Gordon or, 
equivalently, the massive Thirring model, which possess an infinite number of conservation 
laws. While generic models are not integrable, the integrability of their low-energy theories 
implies that their low-frequency, low-T optical conductivity can be strongly influenced by the 
presence of approximately conserved quantities ||,|^. 

Quantitatively, the influence of conserved quantities (CQ) Pn on the conductivity can be 
understood in the following way. Any linear combination of the Pn is conserved: the P„ span 
a vector space. Components of J "perpendicular" to this space will decay. The component 
"parallel" to it, i.e. the projection Jc of J onto this space, is conserved and will give rise to 
an infinite conductivity. What is a physically meaningful definition of "perpendicular"? We 
define an operator B to be perpendicular to A, if (B) = in a (linear-response) experiment 
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where a field conjugate to A is apphed, i.e. A J- B Xab — 0, where xab is the usual 
linear-response susceptibility. Therefore, the susceptibilities are the natural scalar product of 
operators, {A\B) = Txab- This definition is widely used in memory functional theory Q. 

The total weight of the optical conductivity is given by Re(T(a;)dw = 27r(xjj/2). The 
non-decaying fraction Jc of J will induce a Drude weight, 



where (xc)nm — Xp p is the matrix of susceptibilities of conserved P„, and the projector 

onto the space of P„ is written as Yl,nm 1^") i^c^ l'^)nm (^"^l- (0) ^'^'^ extends over 
all the conserved quantities of the system. If the summation includes only a subset of them 
one obtains a lower bound for D, known as the Mazur inequality Q|. The physical arguments 
given above have been put on a more rigorous basis by Suzuki Q who derived Eqn. (|l]) many 
years ago. 

The interest in (|l|) was recently revived when Castella, Zotos and others [|[|^ realized 
that in integrable ID quantum models the Drude weight D is finite at finite T , induced by 
the conserved quantities which prohibit the decay of the current. In a number of ambitious 
numerical studies of (small) quantum systems [Q-^ , the surprising observation of a finite Drude 
weight at r > in systems like the Hubbard model was qualitatively confirmed. Moreover, 
Castella and Zotos ||] found that the mobility of a single particle coupled to Fermions is 
enhanced close to an integrable point. However, a number of fundamental issues remained 
to be answered, partly due to numerical problems and a lack of a theory to investigate the 
consequences of small deviations from integrability. 

The goal of the study reported here is threefold: first, we want to test the important exact 
expression (|l|) for D numerically for a simple system. Previous results have cast serious doubts 
on its validity 0. Our second objective is to investigate how weakly violated conservation laws 
influence the conductivity and to test analytical approximations for a{ijj). A third goal is 
to calculate the optical conductivity close to a Mott transition in a ID system. This is directly 
relevant to recent experiments in Bechgaard salts fictl where well defined low-frequency peaks 
have been observed in (j{uj) below a Mott gap. 

The derivation of an effective model of charge excitations close to a Mott transition starts 
from the Luttinger liquid description of a nearly half-filled band in the presence of Umklapp 
scattering. Using a Luther-Emery transformation and a subsequent refermionization (see 
e.g. ||l^ for details) the charge-sector maps onto a two-band model of Fermions with a band 
gap 2A, i.e. a "semiconductor". The interactions of the "particles" and "holes" are repulsive 
for 1/2 < K < 1 where K parameterizes the strength of interactions in the Luttinger liquid. 
Any finite repulsive interaction in ID leads to a total refiection of particles in the low-density, 
low-energy limit. Therefore, Sachdev and Damle argued (for a similar problem of a gapped 
spin-chain) that the interaction is well characterized by a hard-core model. Exactly at half 
filling, the densitiy of charge excitations is exponentially small for T < A. Their distance is 
therefore exponentially large, much larger than their de Broglie wavelength hj \/2-KmaT and 
for sufficiently small doping and T < A a classical description is possible A simulation 
of the resulting classical models allows to test (0) and to obtain the optical conductivity with 
high precision and small numerical effort. 

We consider a model of two types of particles with the kinetic energy. 




(1) 



all CQs 




(2) 
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Fig. 1 - (a) Current-current correlation function S{t) for T = 0.5, N+/N- = 200/600 and m± = 
1±0.3. In a few collisions S(t) decays down to 2TDint (Q). Due to the approximate conservation of V, 
it takes much longer (inset) to relax down to 2TD (^. Time and distances are measured in units of 
the average particle spacing L/{N+ + N-); in this units the typical time between collisions is ~ \/ m/T 
for low T. A comparison with the analytical approximation and (|l^) (dashed line) confirms our 
interpretation that despite the huge mass difference, the proximity to the integrable point is at the 
origin of the slowly decaying mode. To allow for a better comparison of the numerical data to the 
analytical approximation, we have multiplied all scattering rates in the latter by a factor of « 2 to fit 
the short time scale ((|ll|) is only valid up to prefactors of order 1). (b) Optical conductivity Re o-{ijj) 
(same parameters, doubled system size). Weakly violated conservation laws induce a low frequency 
peak (see inset). The Drude peak 2-kDS{uj) due to momentum conservation (0) is not shown. 



N± is the number of particles with charge ±1 and momentum p±^i. i?kin together with the 
local hard-core repulsion define our model. The particles are distributed randomly on a finite 
ring and the initial momenta are assigned according to the Boltzmann distribution. 

The conductivity is computed from the ensemble average of the current-current correlation 
ftinction S{t—t') — j^{J{t)J{t')) via the fluctuation-dissipation theorem (in the classical Hmit) 

Recr(w) = —S{uj) = 2itD6{uj) + Rea^^Jto) with _D = — lim S{t). (3) 

2T 2T t->oo 

L is the size of the system and the current is given by 

J(i) = gz;+,,(t)-5]z;_,,(i) (4) 

i=l 3=1 

where v±_i — de±{p±^i) / dp±_i is the velocity of the ith particle. Due to conservation laws, see 
below, the dynamics is not ergodic. Therefore, we average over several thousand runs with 
different initial conditions. In each run, each of the 500 to 1800 particles undergoes typically 
1000 collisions. 

In the following, we compare our numerical results for D with the theoretical prediction 
(|l|). We first consider the situation to+ ^ m_, where we expect that only momentum, 
P = J2i^ P+,i + J2i P-,i^ energy H and charges N± are conserved. H and A^^^ do not 
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contribute to as xh.j — Xn±j — due to symmetry and therefore 

^ 1 xIp ^ J_{JPl_ ^ TjANf 1 

2xpp 2TL{PP) 2L N+{p^)^ + N^{p^)_ 

where N = and AA^ = N^—N- are the total number of particles and the total charge, 

respectively. Conveniently, in our model all static susceptibilities are exactly given by their 
value in the non-interacting theory as we are considering only local interactions in a classical 
system, e.g. xjp = tz {^+ {P'^+)+ ~ (P^-)-); where the one-particle expectation value 
(...)^ in the absence of interactions is given by (A)^ — J dpA{p)e~''^^'^P^/^ / J c?pe~^±^P^/"^. 

The situation is qualitatively different when the masses of the "particles" and "holes" are 
equal, m+ = m_. The system is integrable and more conserved quantities exist. Jepsen 
calculated the properties of this system analytically and we used his results fl^,^ for <j{uj) 
to check our numerics in the non-relativistic limit T <^ m. The whole momentum distribution 
of the particles is conserved for m+ = m_ since in any collision of particles with equal masses 
the momentum is just interchanged. Hence, all moments 

i=l j=l 

of the momentum distribution are conserved. What is the projection of J onto the space of 
CQs? Consider the total velocity V = X^i^^ which is conserved for m+ = m_ and 

differs from J by just a minus sign. Since for equal masses xjp„ = ^^tl^^ ('^p")± — T^XvPn; 
the projection Jc of J onto the space of conservation laws P„ is parallel to V and 

2xvv^ ^TL (VV)- 



For m+ = m_, this evaluates to 2Tl^^W^ (^^) 

In Fig. 1^ we compare the analytical formulas (||) and (|^) to our numerical results. Within 
the statistical errors the computer experiments clearly confirm the theory. The conservation 
laws for equal masses to+ ~ m„ are very different from those for m+ ^ m_. Therefore, D 
is not a continuous function of the mass difference D{Am 0) 7^ Z^int = D{Am ~ 0) as is 
shown in Fig. |^a (dashed line). 

What will happen if Am is small? This is a physically relevant situation, e.g. for a Mott 
insulator with a small gap A <C e_F, where corrections which break the particle- hole symmetry 
are small but finite. Due to this small mass difference Am, the higher moments P„ in (|6|) 
with n > 1 will slowly decay and the following physical picture emerges for |Am| <^ m = 
(m+ + m_)/2: we decompose the current J = Jf + Js (as before) into a "fast" component J/ 
perpendicular to the currents Pn and the projection Jg of J onto the space of slow variables 
P„. After a few collisions, i.e. after the time tf, Jf has decayed completely and S{t) = 
(J(<)J(0)) « 2TDint for tf <^ t ts as is shown in Fig. |^a. Js is exactly conserved for 
Am — 0. For finite Am the projection Jc of Js = Jc + Jsd onto momentum P = Pi is 
still conserved but the component Jsd perpendicular to P decays on the scale ts oc l/(Am)^. 
Accordingly, S{t) decays further to S{t ^ ts) = 2TD, where the decay stops (Fig. |l|a). 

The long-time behavior of S(t) determines the low-frequency behavior of the optical con- 
ductivity (Fig. I^b): at zero frequency, Rea(a;) is characterized by a 5-peak with weight 2TrD. 
For small w, we expect a sharp peak of width F^ = 1/ts, weight 27r(Dint — J?) and with a height 
approximately given by (I?int — D)/rs. The slow decay rate determines Re 17(0; ^ 0) M. 
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Fig. 2 - (a) The numerically determined Drude weights D{T) both for an integrable model with 
m+ — rri- = 1, (diamonds) and a non-integrable model (m+ — 1, m_ — 2.7) (triangles) agree 
within the statistical error bars with the theoretical prediction (^) and (^) (thick and thin solid lines) 
(A'^^ = 3N-). The dashed line shows the Drude weight (^) of a non-integrable system in the limit 
I Am I —» 0''', m± ~ 1. (b) Ratio of the fast and slow decay rate Ff/Fs as a function of Am for 
m± = 1± Am/2, N+/N- — 100/300 and T — 0.25. Close to the integrable point, the slow decay rate 
vanishes oc (Am)^. The points are obtained from a fit to S{t) determined from numerical simulations. 
The analytical result ( pl]Jl^ ) (line) describes the observed behavior qualitatively correct. The good 
agreement at high Am is accidentially, ([ll] ) is expected to be valid only for small Am and only up 
to prefactors of order 1. 



In the following, we try to give a crude analytical estimate of the slow decay rate Tg within 
the framework of the (classical) memory matrix formalism [Q . The correlation function 

Cmnit) = (j™(t)j„(0)) (8) 

of a set of currents jQ,ji, ■■■,jN is calculated from 

{M{lo)C-' - iu;)C{lo) = C with M™„(a.) = {dtjm \ Q | 9*J„) (9) 

where Cmn = C{t = 0),nn = Txj^j„ is the matrix of susceptibilities, M the so-called memory 
matrix and L the Liouville operator defined by the Poisson brackets with the Hamiltonian 
iLA = {A, H}. The time evolution of M is determined by the projection of L onto the space 
perpendicular to the currents jk with the projection operator Q = I — J2ki \jk){C^^)ki{ji\- 

The main idea of the memory matrix formalism is that approximations for M are much 
less "dangerous" than approximations for e.g. C{uj), at least if all relevant slow variables are 
included in the space of the j„ • In this case, the dynamics in the perpendicular space defined 
by Q is dominated by fast processes and therefore M is - hopefully - non-singular in the low- 
frequency limit. Here, we test the applicability of this approach in a situation governed by a 
clear separation of time scales. It is easy to check ^ that indeed the exact Drude weight (P 
is reproduced if all relevant conservation laws are included, e.g. jo = J, jn = Pm n — I...00. 
Actually, it is sufficient to keep track of J and its projections onto the P„ and P. For Am — 0, 
the projection Jg of J onto the space of the P„ is parallel to V with 
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For small Am we replace m by {m^ + to_)/2. We can therefore restrict our analysis to the 
three-dimensional space of operators jo — Ji ji = P Sind j2 — V. One can check that this 
approximation is not only valid for | Am| <C ?n at any temperature but also for T <C m ± Am 
for arbitrary Am because Jg is a linear combination of P and V in the latter case. 

While the weights of all features in the optical conductivity can be calculated exactly, it 
is much more difficult to determine the decay rates for small but finite Am. To estimate the 
various elements of the memory matrix, we use two approximations: we neglect all correlations 
between collisions and we neglect the projection operator Q in (||): 

M„™(t) « (a,j„(t)%„(0)>«^(AjWAjW,5(t,))<5(t) 

i 

^ ^^(Aj„Aj„>+-«_|)p^p_<5(t) (11) 

where Aji'-* is the change of j„ at the ith collision at time ti, dtjn — J2i ^i'n' ^{t — ti) 
and(/p+,p_)p p_ =//p,+p_e-(<^+(f+)+^-(P-»/'^//e-(^+(P+)+^-(P-»/'^ is the thermal average 
over the momenta p± of the two incoming particles with velocities v± and charge + and 
respectively. Note that collisions of particles with equal charge do not change j„. Both 
approximations in are uncontrolled and will induce errors of order 1 in the prefactors 
of the collision rates. Nevertheless, a comparison with our numerical data shows that the 
qualitative dependence of cr as a function of T, Am and AN/N is very well described by ([lT|). 
For = N- and Am — and T ^ we obtain e.g. cr{T) = ^/tt / {l&^J mT) to be compared 
to the exact result ||l2,13| (t(T) = l/-\/27rmT, the error is smaller at higher T. 



The size of Aj„ determines the various collision rates. Ajo ~ 2(w(p+) — v{p-)) is finite 
for Am and Mqo describes the initial fast decay of J. Aji vanishes due to momentum 
conservation. Accordingly, the projection of J onto P will not decay. Aj2 is linear in Am as 
V is conserved for Am = which leads to a slow decay rate proportional to (Am)^ 

r./F/ oc (Am)2. (12) 

A comparison of the analytical result with the numerical data shows that our analytical calcu- 
lation within the memory matrix formalism does not only reproduce the exact weights (Fig. |^) 
but also describes the relevant decay rates Fg and Fj qualitatively and semi-quantitatively 
correctly. In Fig. |^b, we compare the exact decay rates determined from a fit to the numer- 
ically calculated S{t) ~ cie~'"^* -f C2e~'"=* -I- C3 to the corresponding rates derived from our 



analytical approximations (H) and (11). In the limit Am ^ we obtain 



M22 - M^^/Moo ^ M( 



00 



C22 - Cf,/Ci, ' ' ^ Coo(l - (A7V/iV)2) 

where the memory functions M„,„ are evaluated at w = 0. At low T we get e.g. Tf fa S^/T/tt 
and r, « 6(Am)2(l - {AN /Nf)y^Tj^. 

Within our (crude) approximation scheme, S{t) decays exponentially to a constant for long 
times. This is certainly not correct, as e.g. for Am — 0, S{t) ~ const. — 1/t" with a = 3 
as is known from the exact result [|l^,|l^. Our numerical results suggest that the (algebraic) 
long-time tails are more pronounced for Am 7^ but we were not able to extract the exponent 
a reliably due to large finite size effects (due to the hard-core interaction the particles cannot 
pass each other and their diffusion is stopped after a distance of order \/]V). The main physics 
discussed in this paper, i.e. the existence of a time-scale l/F^ cx l/(Am)^, seems not to be 
affected by the presence of algebraic long-time tails. 
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How small has Am to be that (j{u!) displays a well defined peak of width Tg due to the 
slowly decaying component of J? To decide this question, we compare Re (7(0; Am 
0+) w 2ADs/Ts with Recr(w ^ 0+,Am = 0) « 2ADf/Tf where AZ?, = Ant - D and 
AZ)/ = XJjf^ ~ -Dint are the weights of the slowly and fast decaying modes, respectively. The 
low frequency peak is visible if 

ADf ,AD, ^ {Amfjl [AN/Nff 

' V, ^ 2T^{AN/NY ^ ' 

where the last equation is valid for low T. The peak is not observable for T < |Am| which 
is a consequence of the fact that the model is Galileian invariant for T — > 0. The slowly 
decaying modes are much less relevant for AA^ — 0, as xpj — AN/L and no well pronounced 
low-frequency peak develops for low T j^]. 

In this paper we have investigated numerically and analytically the optical conductivity 
for a model of charge excitations above a Mott gap. Due to conservation laws, the Drude 
weight is finite at T > for AA^ ^ 0. We can confirm predictions (0,|^) for the Drude weight 
and (T(tj) both for an intcgrable model and in a situation where conservation laws are weakly 
broken and low-frequency peaks arise in cr{uj). One of the surprises of this study was that 
even for a really large mass difference, e.g. m± = 1 ± 0.3, the proximity of the Hamiltonian to 
some integrable point with m+ = m_ strongly influences the low-frequency conductivity (see 
Fig. |l|b)). At least for the model studied here, the memory matrix formalism is a reliable tool 
to study this type of physics qualitatively. Furthermore, we expect a quantitative agreement 
in situations where it is possible to calculate the memory matrix perturbatively. One may 
speculate that the low-frequency peaks in the optical conductivity of the Bechgaard salts are 
related to approximate conservation laws as they are characterized by tiny weights, unusual 
line-shapes and a very long mean- free path . 

We would hke to thank N. Andrei, K. Damle, F. Evers, P. Wolfle and X. Zotos for helpful 
discussions and the Emmy-Noether program of the DFG for financial support. 
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